;$Id: lognormal.pro 6 2009-07-01 08:17:54Z dingyy.pku $
FUNCTION LogNormal, N, d_bar, d, sigma
  pi=3.14159265D;
  logSigma=Alog(sigma)
  twoPiSqrt=Sqrt(2.0d*pi);
  result=N/(twoPiSqrt*logSigma)*Exp(-(Alog(d)-Alog(d_bar))^2/(2.0d*logSigma^2));
  Return, result;
END
